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ON THE STABILITY OF THE MOTION OF A VISCOUS LIQUID* 

BY 

FRANCIS ROBERT SHARPE 

§1. Introduction. Reynolds f has established an energy criterion for the 
stability of the motion of a viscous liquid, and has applied it to the case of a 
liquid moving between two parallel planes. For a liquid of density p and vis- 
cosity /t, moving between two planes 26 apart with mean velocity U he found 
that the motion is unstable when 

?^>517. 

Using the principle that the critical velocity is inversely proportional to the 
hydraulic mean depth he inferred that for a cylindrical pipe of radius a the 
critical velocity is given by 

^^>1034, 

and he compared this with his experimental value 1900. In the following paper 
I discuss directly the case of the cylindrical pipe and find that the motion is 
unstable when 

2apU^ 



fJ' 



>470. 



I also find, by using a different solution of the equation of continuity and the 
boundary conditions, that the motion between two parallel planes is unstable when 

instead of 517 as found by Reynolds. 

§ 2. JRectilinear flow through a pipe. When the motion through a pipe is 
symmetrical the equations of motion, expressed in cylindrical coordinates r,d,z, 
are J 



* Presented to the Society February 25, 1905. Received for publication May 5, 1905. 
fREYNOLDS, Philosophical Transactions of the Royal Society, vol. 186 (1895) 
A, pp. 123-164. 

J Love, Mathematical Theory of Elasticity, vol. 1, p. 217. 
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du du\ dp^ dp^ Prr-P,9 



(1) p|^_+„_ +„_j=_ +-5r + 



r 



dw dw dw\ dp^ dp p^ 



,„. fdw dw dw\ d 



5r 5z r * 

The extensions are* du/dr, u/r, dw/dz. and the shear du/dz+ dw/dr. 
Hence f 

( du dw \ 
The equation of continuity is J 

/OS ^^ « ^^ /v 

(3) ar + -r + -^ = <^' 

and the boundary conditions are 

(4) M = w = when r = a. 

The simplest solution of these equations is § 

dw dw dp dp 

u == ^TT- = -^— = -5- = , -5- = constant, 
c)< 5z 5r 5a 

- = -4^S(«^-'^)' 
and the mean velocity of flow is 

8/t dz ' 

§ 3. Steady non-reetUinear motion. For some minimum value of CT", to be 
determined, permanent eddying motion is possible, that is, the rectilinear motion 
is unstable. The mean motion is still rectilinear but u is no longer zero at any 
particular instant. We assume then that 

(5) u = u', w = w+w', 

where u) is a function of r, and u' and w' are periodic functions of r and z whose 
mean values are zero. 

The method used by Eetnolds is to calculate the average rate of increase of 
the motion m = u', m> = to', which is the motion relative to the mean motion 
M = , w = w, for a volume of liquid extending to the boundaries of the liquid. 

*Lovk, loc. cit., p. 216. 

t Lamb, Hydrodynamics, p. 512, 521. 

% Lamb, loo. oU., p. 134. 

i Lamb, loo. oit., p. 521. 
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§ 4. 7^ equations of mean motion. Using a bar to denote a mean value 
we have, on substituting (5) in (1), (2) and (3), 

) 

P„+P'rr—P»»—P'»» 






+ 



{dw ,dw' ,_ ,x ^w' ,5io\ ^ip„ + p'^ 

, ^iPz,+P'..) , Pr,+P 'r. 

+ dz + r ' 

du' u dw' 

<6) w+v + w-^' 

Multiplying (6) by pu' and pw' respectively and adding to the preceding equations, 
(dw dju'u') d(u'w') du' u'u' I a(^^+;>:j 

, ^{P,r + P',r) , P»-^P'rr-P »i, - P'»e 

\dw' d{u'w') d{w'w') dw' dw u'w 1 5(p„-KKJ 

■^ az + r • 

Taking mean values, we have the equations of mean motion 



/7\ . d(u'u') U'U'\ 5Prr,^P.r,Pr, 



/«\ „l^(«V), «'«''! ^Pr.,^P^,Pr. dp fd'^w ldw\ 

^^^ ''\~W~ + ~r~ | = 'ar"*"ar + 7 = ~a^ + ^V^+ra;;- 

Multiplying (8) by r and integrating 

— — . r" a» aio 

Since this holds for r = , we have C = . Hence 

dw r dp -,—f 

§6. 2%e energy of the relative motion. Subtracting (7) and (8) from the 
preceding equations we have for the relative motion, 
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{du' ^ , , , —r—r^ _ 5m' d(u'w'\ u'u' —u'u \ 

_^JPrr,^jL . P'rr-P'ti, 

~ dr '^ dz '^ r ' 
idw' ,dw' ^ , , , — — 7v _ 5m)' dw ■' 

P^r + "^;r + ^(««' -uw) + w-^ + u -^ + 



dt ^ dr ^ dr^ — /^"^ dz ^ dr 



^P'r, ,^P'» , P'r. 

dr dz r 



Multiplying these equations by u' , and w' respectively, adding, and taking mean 
values, 

valueof.'(^^^^^^^)^.'(^=^^^^). 



mean 



Integrating over the section of the pipe and by parts on the right hand side, 
dividing by 27r, and remembering that «', w' vanish on the boundary, we 
obtain the equation 

J« [dt'^'^dz)^^^'' ■•"'* )rdr = J^-puw -^rdr 

where, in the second integral on the right, mean values are to be taken. 

The left hand side of this equation is the average rate of increase of the energy 
of the relative motion of a length l/27r of liquid in the pipe. 

Hence we have Reynold's criterion for the possibility of permanent eddying 
motion expressed in cylindrical coordinates, namely, 

(10) -] ,„'.' -^rdr^ ,[ [2 (-,-)+ 2 (^ ) 

u'^ (du' dw'y\ , 

where, in the integrand on the right, mean values are to be taken. 

The interpretation of equations (9) and (10) is that in order to maintain the 
relative as well as the mean motion an increased pressure gradient is necessary 
and the additional energy is dissipated in the relative motion, the right-hand 

Trans. Am. Math. Soc. 34 
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integrand of (10) being Stokes's dissipation function expressed in cylindrical 
coordinates. 

§ 6. Determination of a possible relative motion. The solution of the con- 
tinuity equation (6) with the boundary conditions u' = w' = 0, for r = a, 
which leads to a minimum value of U is 

Try *jrlj% trv itIjz 

u' = LA^^ (sin p + sin 3jo) sin -^ — LB ^^ (sin 2p + J sin ip) cos -^ , 

2(sin2j + sin 3/>) + „- (cos^) + 3 cos3/)) cos-^— 

{trr 
2(sin 2/) + J sin 4p) + 3- (2 cos 2/> + 2 cos 4p) 



. irLz 
^^°-2j' 



as is easily seen on substitution, p being irrj^a . 
These give 

LABir^r^ 

u'w = — qo~2 — (^ sin^ + 3 sin %p — sin hp — sin 7/)). 

Hence from (9) we find, neglecting m'w'*, the left-hand side of equation (10) 
equal to 

LABir'pU 



64a= 



'pU r 

^-— I r<(3sin^-|-3sin3/) — sin5j9-sin7j9)«Z»- = .216Z^5apC;'. 



Substituting the values of u' and w' in the right hand side of (10), taking mean 
values, and integrating, we obtain the condition 

•^1^^:^^^^(.304^^ + .1275^)Z^ 

+ (4.26^^ + 6.675^)i>^ + 35.7^2 + 16.34£^ 
Hence for a minimum value of ZZwe have to minimize 



(■f..mx)^..(i^%e..,.)^.(?f'.ia.M.)-i 

where \ = 5/ J. . 

Differentiating with respect to L and \ we find 

3 (.304 + .127X')i^ + (4.26 + 6.67X*)i' - 36.7 - 16.34\* = 0, 
(_ .304 + .127X^)Z^ + (- 4.26 + 6.67X')i^ - 35.7 + 16.34\2 = 0. 
Eliminating X' , 

.4633Z' + 20.55i/« + 151.67Z* - 2333.35 = 0. 
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Hence 2/ = 1.798, \ = 1.158, approximately. The minimum value of U ia 
therefore given by 

-^^^ = 470. 

Reynolds' experiments gave 

Reynolds' result for the flow between two parallel planes a distance 2b apart is 

A* 

and because the hydraulic mean depth for a cylinder radius a is |a and for two 
parallel planes is 6 he inferred that for a cylinder 

?-°e?_ ,034. 

§ 7. Flow between parallel planes. When the motion is the same in planes 
parallel to the xy plane the equations of motion are * 

/a« a« atj\ dp, dp^^ 



where 



(du dv\ 
% + aij' 

The equation of continuity is 



du dv 

dx dy 

In the case of liquid flowing between two parallel planes y = zkib the bound- 
ary conditions are u = v = when y = ±b. 
The simplest solution of these equations is 

du dp dp 1 dp 

and the mean velocity of flow is 

3/i dx 



* Lamii, Hydrodynamics, pp. 512-520. 
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Using the method of the preceding problem we find for the equation of mean 
motion 

d(u'v') dp d^u 

(11) ''^--i + ^F^-' 

and the criterion of stability, 

§ 8. Determination of a possible relative motion. The solution of the con- 
tinuity equation, with the boundary conditions, which leads to a minimum value 
of Uia 

u = A(^Aap 4- sin Sp) cos-^r — I- B{2, sin 2p + 4 sin 4p) sin „, , 

v' = — LA(cosp + ^cos3p)sin-2v- + LB (cos 2p -f cos4p)cos-^Y- » 

where ^) = iryl'ih . 
These give 

4m' ■«' = — LAB '{ Jg*- sin 2) + 6 sin 3p + | sin 5p + \ sin Ip } . 

To determine du^ dy we have, on integrating (11), 

—r-, dp du 

Now u'v' is an odd function of y and so is du/dy since the mean motion is 
symmetrical with respect to the planes y = zhb. Hence C is zero and approxi- 
mately 

du y dp SUy 

dy~ fi dx~ 6* ' 

Hence the left hand side of (12) reduces to 

- 18.6 LABpU 

Substituting the values of u' and v' in the right hand side of (12), taking mean 
values, and integrating, we have the condition 

-IS.GLABpU 

-^{'if(-V-^'+2^') +^'(2^^+205') + ^z(10^H2725») }, 



or 



BbpU^ /^TT V(-V-7l' + 2B'-)L' + (4^' + ^(iB')U + 10^' + 272^' 



(i) 



IM ~\2l -3.1LAB 
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Reynolds takes the solution 

u =^(cos^ + 3 cos 3/))cos-o^ + -B(2 cos 2^ + 2 cos 4p) sin g, , 

v' = LA(^sinp + sin 3^?) sin-gr LB {sin 2/> + J sin 4p)co3 q, , 

and obtains the condition 



^Ki)' 



UUp^(nr\ (2^' + f^')X^ + (20^' + 16^^ )i:' + 82^.' + 80^' 

-1.3252/^5 • 



The coefficients of A^ and B'^ being nearly equal, he determines the minimum 
value of £7" by assuming A = B. 

It is interesting to notice that his value for v' vanishes midway between the 
planes and u is there a maximum whereas the values here adopted give a maxi- 
mum value of v midway between the planes, which seem to be more in accord- 
ance with experiments. 

Using the method of the preceding problem we find a minimum value of U 
given by 

when L = 1.398 and AjB = 4.03. 

Hence our final result is that the motion is unstable when 

-gl67. 

M 

Reynolds' result is 517. 

COBNKLI, UNIVKRSITY, April, 1905. 



